Direct evidence of overdamped Peierls-coupled modes 
in TTF-CA temperature-induced phase transition 
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In this paper we elucidate the optical response resulting from the interplay of charge distribution 
(ionicity) and Peierls instability (dimerization) in the neutral-ionic, ferroelectric phase transition of 
tetrathiafulvalene-chloranil (TTF-CA), a mixed-stack quasi-one-dimensional charge-transfer crystal. 
We present far-infrared reflectivity measurements down to 5 cm -1 as a function of temperature above 
the phase transition (300 - 82 K). The coupling between electrons and lattice phonons in the pre- 
transitional regime is analyzed on the basis of phonon eigenvectors and polarizability calculations 
of the one-dimensional Peierls-Hubbard model. We find a multi-phonon Peierls coupling, but on 
approaching the transition the spectral weight and the coupling shift progressively towards the 
phonons at lower frequencies, resulting in a soft-mode behavior only for the lowest frequency phonon 
near the transition temperature. Moreover, in the proximity of the phase transition, the lowest- 
frequency phonon becomes overdamped, due to anharmonicity induced by its coupling to electrons. 
The implications of these findings for the neutral-ionic transition mechanism is shortly discussed. 



I. INTRODUCTION 

In a formulation similar to a student's exercise, more 
than fifty year ago Rudolf Peierls pointed out that one- 
dimensional (ID) metals are intrinsically unstable with 
respect to lattice distortions opening a gap at the Fermi 
energy. 1 It did not take long to realize that the quasi- ID 
electronic structure of organic charge-transfer (CT) con- 
ductors offers an almost ideal model for the Peierls mech- 
anism. Extensive investigations in the 1970s revealed 
that the Peierls instability may occur also in ID Mott 
insulators, due to coupling of the lattice to spin (spin- 
Peierls), 2 or to both electronic and spin degrees of free- 
dom ("generalized" Peierls instability). 3 At the time, the 
attention was almost exclusively devoted to segregated- 
stack CT crystals, where identical electron-donor (D) or 
acceptor (A) 7r-molecules are arranged face-to-face along 
one direction. Mixed stacking (...DADA...) received far 
less attention, at least until a new kind of phase transition 
was discovered in these systems: the so-called neutral- 
ionic phase transition (NIT) It is characterized by a 
change in the ionicity g, i.e. the charge transfer from D to 
A in the ground state. Temperature T or pressure p may 
induce a change from a quasi- neutral TV, g < 0.5) to a 
quasi-ionic (/, g > 0.5) ground state. It was early recog- 
nized that the ionic state is subject to the Peierls insta- 
bility, yielding a stack dimerization (..DA DA DA..)P^ 

NIT are complex phenomena, resulting from the inter- 
play between two order parameters, the ionicity g and 
the dimerization and present a rich and intriguing 
phenomenologyP^ Some of the phenomena occurring on 
approaching the phase transition, like the dramatic in- 
crease of the dielectric constant or the observation of 



diffuse scattering signals in x-ray 

diffractionJ^EU have 
been recently ascribed to charge oscillations associated 
with the soft mode that yields the stack dimerization, 
namely, the Peierls modeP^ 13 * 14 * It is then of fundamen- 
tal importance to identify the Peierls mode in order to 
confirm the above suggestion and to shed light on the 
NIT mechanism. At variance with segregated stack crys- 
tals, the Peierls mode is optically active in mixed stack 
structures^ which offers a rather unique possibility for 
its optical detection and characterization. Finally, since 
the dimerized ionic stack is potentially ferroelectric,^ the 
Peierls mode is the counterpart of the so-called ferroelec- 
tric mode of perovskite-type crystals.^ 

Direct and clear experimental identification and char- 
acterization of the Peierls mode and of its role in the NIT 
has proved to be more difficult than expected! 17 * 181 One 
of the reasons is that the phonon structure of molecular 
crystals is very complex, and several modes are likely cou- 
pled to CT electrons. Therefore the analysis of the spec- 
tra requires a detailed and reliable understanding of the 
lattice-phonon dynamics. We have recently undertaken 
such analysis for tetrathiafulvalene-chloranil (TTF-CA), 
which undergoes a first-order temperature induced NIT 
at 81 KP In a first paper (hereafter Paper I)P^ we pre- 
sented a quasi harmonic lattice dynamics (QHLD) calcu- 
lations of TTF-CA lattice phonons above and below the 
transition temperature, and gave a first partial interpre- 
tation of the infrared (IR) and Raman spectra. 

The specific goal and achievement of the present paper 
is the identification of the Peierls (soft) modes in TTF- 
CA neutral phase. We present detailed and improved 
data for single crystals of TTF-CA in the spectral region 
5 — 200 cm -1 collected in the temperature range down to 
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82 K. Particular attention is devoted to the submillimeter 
spectral range, below 30 cm -1 . The experimental data 
are analyzed in terms of a multi-mode Peierls model, that 
combines the information on lattice phonon dynamics re- 
ported in Paper I, 19 with that on the electronic struc- 
ture as obtained from the exact diagonalization of the 
ID Peierls-Hubbard Hamiltonian in Ref. [13j The model 
fully accounts for the experimental data and allows us 
to demonstrate that mode mixing and overdamping are 
the key to explain the changes of the far-IR spectra on 
approaching the transition. The effective Peierls mode 
detected by an analysis of the combination modes in the 
mid-IR region 20 is also accounted for. In addition, the 
dielectric constant anomaly at the transition^ is quanti- 
tatively reproduced. 



II. EXPERIMENTAL METHODS 

TTF-CA single crystals were grown by sublimation. 
Crystals as large as 4 mm x 4 mm allowed for measure- 
ment down to 5 cm -1 . The reflectivity was probed on 
well-developed naturally grown ab crystal faces along the 
optical axes. The crystals were oriented at room tem- 
perature according to the known spectra in mid-IR re- 
gion, using a Hyperion IR microscope attached to Bruker 
IFS66v Fourier-transform spectrometer. In this paper we 
present the spectra polarized parallel to the a stack axis; 
only these spectra are relevant to the analysis of Peierls 
coupling in the neutral phase 

The reflectivity spectra were collected by three dif- 
ferent setups and techniques, depending on the spec- 
tral region of interest: (i) In the 5 — 30 cm -1 region, 
we employed a quasi-optical submillimeter spectrome- 
ter, equipped with backward wave oscillator as coherent 
and tunable source; a Golay-cell served as a detector. 
Temperature dependent measurements were conducted 
in an exchange-gas helium cryostat. The absolute reflec- 
tivity was obtained by comparing sample reflectance to 
reflectance of an aluminum mirror fixed on the same aper- 
ture. The spectral resolution in this range is better than 
0.01 cm -1 , (ii) In the 20 — 150 cm -1 range the reflec- 
tivity was measured by a Bruker IFS66v with a spectral 
resolution of 0.5 cm -1 . The spectrometer was equipped 
with 1.4 K Si bolometer and a CryoVac cold- finger cryo- 
stat. The sample was fixed on the cold finger by car- 
bon paste; a good temperature contact was ensured by 
Apiezon paste. To get absolute values of reflectivity, the 
in-situ gold-evaporation technique was adopted: 21 after 
measuring the temperature dependence of the sample re- 
flectivity, gold is evaporated onto the sample surface, and 
used as reference, (iii) Reflectance measurements in the 
100 — 6000 cm -1 range were done in a Bruker IFS 113v 
FTIR spectrometer equipped with a 4.2 K Si bolometer, 
an MCT detector, and an exchange-gas helium cryostat. 
The absolute reflectivity values of the sample have been 
obtained by a comparison to a Au mirror of the same 
size. In this spectral range a resolution of 2 cm -1 was 



chosen. The reflectivity values agree well with available 
literature data, 22 which extend from 3000 to 14 000 cm -1 
and thus cover the charge transfer transition. The exten- 
sion to high wavenumbers is needed to obtain reliable 
Kramers-Kronig transformation (KKT) of the data. 

We ensured that good overlap exists between the data 
measured in the three spectral regions. In order to match 
the absolute values, the submillimeter data were rescaled 
to those obtained by the Bruker IFS66v FTIR because 
the in-situ gold-evaporation technique ensures the most 
precise measurements of the absolute reflectivity and 
does not depend on the quality of the sample surface. 

In the KKT analysis the data were extrapolated below 
5 cm -1 to fit the known values of the static dielectric con- 
stant eiP^We also verified that an extrapolation to a con- 
stant value does not significantly affect the portion of the 
spectra for which experimental data exist. At frequen- 
cies above 14 000 cm -1 an extrapolation R(uj) oc uo~ 2 has 
been used. 



III. REFLECTIVITY SPECTRA 



A. Spectral data 
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FIG. 1: (color online) Experimental and calculated (Section 
|VJ reflectivity spectra of TTF-CA at T = 300 K, with the 
polarization E \\ a. The spectral region above 3000 cm -1 
(dashed curve) is taken from Ref. 1221 Notice the logarithmic 
frequency scale. 

Fig. [T] shows the room temperature TTF-CA reflec- 
tivity along the stack axis for the whole measured spec- 
tral range, with extension to 14 000 cm -1 taken from 
Ref. [22j The reflectivity is rather low, as expected for a 
non-metallic compound. The charge transfer electronic 
transition occurs around 5000 cm -1 and at lower fre- 
quencies reflectivity gets down to about 0.2, with no 
relevant feature observed until we encounter the highest 
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frequency out-of-plane intramolecular modes, at 701 and 
682 cm -1 , corresponding to CA and TTF, respectively.^ 
Below 250 cm -1 we find a couple of other modes of pre- 
vailing out-of-plane, intramolecular character, and finally 
the intermolecular modes we want to discuss in detail in 
the present paper. 
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FIG. 2: (color online) Reflectivity spectra of TTF-CA single 
crystal along the a axis taken at different temperatures T > T c 
above the N - I transition. 

The temperature dependence of the low- frequency (5 — 
120 cm -1 ) reflectivity spectra is presented in Fig. At 
room temperature a group of strong phonon bands is ob- 
served just below 100 cm -1 . At lower frequencies the re- 
flectivity remains approximately at the level of 0.4 down 
to the lowest measured data (5 cm -1 ). On cooling the 
sample from 300 to 82 K, the reflectivity above 120 cm -1 
reveals only unimportant changes, and is therefore not re- 
ported. Instead Fig.[2]shows that dramatic changes occur 
below 100 cm -1 , where reflectivity rises by 1.5 times on 
lowering temperature to 82 K, just before the NIT. As 
demonstrated in Paper Ip^ below the transition temper- 
ature (T c = 81 K) the reflectivity drops down to values 
below those at T = 300 K, and does not change consid- 
erably on further cooling to 10 K. 

A more detailed inspection of Fig. [2] provides evidence 
that the strongest increase of reflectivity, from about 0.4 
to 0.6, occurs on cooling between 300 and 150 K. Some in- 
crease of reflectivity below approximately 40 cm -1 occurs 
on further cooling, indicating a shift of spectral weight 
to lower frequencies. The quality of the data in the low- 
frequency region, essential to appreciate the shift of spec- 
tral weight, is illustrated in Fig. [3j where the reflectivity 
at T — 82, 90 and 100 K is normalized to the reflectivity 
at 110 K. This figure illustrate the good overlap between 
submillimeter and FIR data, collected with two different 
instrumentations (see Section |n|) . The reflectivity ratio 
R(T)/R(110 K) for these two spectral regions have in- 
deed differences smaller than the noise of the measure- 
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FIG. 3: (color online) Ratio of the of the reflectivity spectra at 
82, 90, and 100 K to the reflectivity at 110 K. Measurements 
in the range below 25 cm -1 performed by using the submil- 
limeter spectrometer (points), data above 25 cm -1 collected 
by using Bruker 66v spectrometer (lines). Note the increase 
of reflectance below 25 cm -1 at temperatures close to the N 
- I transition temperature. 



ments. Thus we safely confirm that by lowering T below 
100 K the reflectivity drops in the range 30-100 cm -1 
range and increases below, with a redistribution of the 
overall spectral weight. 



B. Optical conductivity 

Frequency dependent conductivity spectra, obtained 
by the KKT of the reflectivity, are presented in Fig. [4] 
In order to better understand the temperature evolution 
of the conductivity spectra, and to gain confidence in 
the band assignment, we have also fitted the original re- 
flectivity data with a set of Lorentzians. The resulting 
conductivity practically coincides with the KKT conduc- 
tivity, as exemplified for the 300 and 82 K data in the 
upper panel of Fig. [5] The individual Lorentzians used 
in the deconvolution are also shown. In the bottom panel 
of the Figure we report an enlarged view of the temper- 
ature evolution of the KKT conductivity vs. the multi- 
Lorentzian reflectance fit. The logarithmic scale for the 
wavenumber axis focus attention on the deconvolution of 
the lowest frequency phonon mode, that we shall discuss 
in detail below. 

Six peaks are clearly identified in the spectral region 
below 100 cm -1 . We label them v\ to v& in order of in- 
creasing frequency, as shown in the upper panel of Fig. [5j 
The temperature evolution of the v\ to v§ peak frequen- 
cies is reported in Fig. [6] In a group of bands with a broad 
maximum around 68 cm -1 at T = 300 K, we resolve three 
modes (^6,^5 and v^) at lower temperatures. While the 
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FIG. 4: (color online) Conductivity spectra along the a stack 
axis at temperatures above the N - I transition. An error 
bar for the absolute values of low frequency conductivity is 
shown. 



high-frequency band shows normal hardening on lowering 
T, the ^5 and v± modes essentially do not change their 
positions. The v<$ band at 47 cm -1 exhibits a very weak 
softening of 2 cm -1 below T = 130 K. The deconvolu- 
tion of the low frequency region of the room temperature 
spectrum suggests the presence of bands around 32 and 
17 cm -1 (y<i and v\), "appended" to the 47 cm -1 band. 
When going from ambient T to 200 K, the v\ band shifts 
below 15 cm -1 and becomes wider. It shifts to even lower 
frequencies and broadens as temperature is reduced, get- 
ting down to about 5 cm -1 for T = 90 K (Fig. [6j, as- 
suming the shape of a broad background at temperatures 



bottom panel), 
show a huge in- 



just above the phase transition (Fig. [5 

The conductivity spectra in Fig. 4 
crease of overall intensity by lowering T. A quantitative 
estimate of this growth in terms of the spectral weight 
between and 100 cm -1 is shown in Fig. [7fi. The inten- 
sities of the individual bands, estimated by integrating 
in the relevant spectral region, are reported in the lower 
frame. Fig. ^f)p shows that only some bands gain inten- 
sity by lowering T: the and v± (estimated as a whole 
maximum around 70 cm -1 since the bands are not well 
resolved), the v<$ band at 47 cm -1 , and the v\ band. The 
intensity of the latter increases only below 100 K. 

While the spectral weight increases on cooling, it also 
redistributes towards lower frequency bands. The tem- 
perature dependence of intensity around 70 cm -1 has 
smaller slope below 200 K, whereas the spectral weight 
of us (47 cm -1 ) band saturates at temperatures below 
130 K. At this temperature the intensity of the v\ band 
starts to rise rapidly. Just above the transition, from 90 
to 82 K, the v\ band grows at the expenses of the v% one. 

With lowering temperature the v<$ and v\ bands show 
considerable anharmonicity. While at T = 300 K the v<$ 
band (47 cm -1 ) can still be fitted with a Lorentzian, it 
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FIG. 5: (color online) Upper panel: Comparison of the TTF- 
CA KKT conductivity at 300 K (green line) and 82 K (orange 
line) with the fit of reflectivity spectrum by a minimum num- 
ber of harmonic oscillators. The Lorentzian bands used in 
the deconvolution are shown as dashed lines, and the result- 
ing conductivity as blue lines. Bottom panel: Comparison of 
the KKT conductivity with the multi-Lorentzian fit at two 
different temperatures. The Lorentz component of the lowest 
frequency phonon (y\) is also shown. Notice the wavenumber 
logarithmic scale. 



becomes wider and very asymmetric at temperatures be- 
low 250 K, as its lower frequency wing grows on cooling. 
The temperature dependence of the full width at half 
maximum (FWHM) of the v<$ band, directly estimated 
from experimental spectra (Fig. [5j bottom panel), is dis- 
played in Fig. [8] The same plot shows also the dramatic 
increase of the FWHM of the lowest frequency v\ band 
below T=100 K, with clear overdamping at temperatures 
close to the transition. 

The present analysis of the spectra evidences a very 
specific behavior of ^5, ^4, and especially ^3 and v\ 
bands as the temperature approaches the NIT. The bands 
clearly correspond to the most strongly coupled Peierls 
modes. By lowering T, they show a huge intensity in- 
crease, with a concomitant intensity redistribution, and 
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FIG. 6: Temperature dependence of the frequencies of the six 
phonons bands observed below 100 cm -1 , labeled v\ to vq 
going from the lowest to the highest frequency. 



important deviations from Lorentzian bandshape. Below 
100 K, we also observe appreciable softening and large 
broadening for the lowest frequency mode. A quantita- 
tive explanation of this complex temperature behavior of 
the spectrum is offered by the model discussed in Sec- 
tion |V| while in the next Section we describe the compu- 
tational methods. 



IV. COMPUTATIONAL METHODS 

The separation of intramolecular vibrations from in- 
termolecular, or lattice, phonons is a common and use- 
ful approximation in dealing with the complex phonon 
spectra of molecular crystals. Thus lattice phonons de- 
scribe translations and rotations of the rigid molecules 
(rigid molecule approximation, RMA). In the framework 
of a Hubbard-model description of the electronic struc- 
ture (see below), the electron-phonon coupling can also 
be separated into two contributions: The totally symmet- 
ric molecular vibrations are assumed to couple with elec- 
trons through modulation of on-site energies (electron- 
molecular vibrational, or Holstein, coupling). On the 
other hand, lattice phonons and possibly out-of-plane 
molecular vibrations are expected to modulate the inter- 
molecular CT integral (Peierls coupling) P3 

In order to characterize the Peierls coupling in TTF- 
CA, we have first calculated the lattice phonon frequen- 
cies and normal coordinates at T = 300 K (TV phase) by 
the quasi-harmonic lattice dynamic (QHLD) method, as 
described in paper iP^The atom-atom potential adopted 
in QHLD accounts for van der Waals and Coulomb forces, 
but does not include the CT interaction. It should 
be noted that in the calculation we have relaxed the 
RMA, by actually considering all the phonons below 
- 250 cm- 1 . 
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FIG. 7: (color online) (a) Temperature dependence of the 
spectral weight of the — 100 cm -1 spectral region (squares), 
compared with the temperature dependence of (OP/dS) 2 
(open stars) (Section [vf. (b) Temperature dependence of the 
spectral weight of the group of bands around 70 cm -1 (1/4 and 
1/5) (rhombs), and of the z/3 (triangles) and 1/1 (circles) bands. 
The lines are a guidance to the eye. 
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FIG. 8: (color online) Temperature dependence of the full 
width at half maximum (FWHM) of the z/i (circles) and of 
the za> (triangles) bands in TTF-CA. 
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We describe the strength of Peierls coupling in terms 
of the linear Peierls-coupling constants: 



9% 



( <%da\ 
2uJi \dQiJ ( 



(1) 



where £da is the CT integral between adjacent DA 
molecules along the stack, and Qi is the normal coordi- 
nate for the i-th phonon with frequency uoi and wavevec- 
tor q = 0. The CT integral and its variation with Qi has 
been calculated by the extended Hiickel method, adopt- 
ing the Wolfsberg and Helmholtz approximation!^ At 
the room-temperature equilibrium geometry, we find t 
= 0.20 eV, in agreement with current estimates.^ The 
gi values, obtained by numerical differentiation, offer a 
reliable indication of the relative magnitude of the cou- 
pling constants. The Peierls coupling strength of the z-th 
phonon is given by Si = (gf/ui), and the total coupling 
strength, or lattice relaxation energy, is Ed = ^Zi^i- 

The electronic structure of TTF-CA is described in 
terms of a modified Hubbard model with adiabatic cou- 
pling to molecular and lattice vibrations. Real space di- 
agonalization of the modified Hubbard Hamiltonian rele- 
vant to stacks of up to 18 molecules with periodic bound- 
ary conditions yields reliable information on the ground 
state properties.^ The expectation value of the dipole 
moment operator and its derivatives are calculated by 
the Berry-phase approach.^ 



V. SPECTRAL MODELING 



As discussed in Paper and elsewhere, 1 171261 in the 
N phase of TTF-CA the Peierls-coupled phonons trans- 
form as the B u species in the C\ h crystal symmetry, and 
are IR active with polarization along the stack axis. We 
shall therefore restrict our attention to these phonons. In 
Table [I] we list the QHLD calculated B u frequencies for 
the equilibrium structure at 300 Kp^and the correspond- 
ing coupling constants and coupling strength obtained as 



described in Section IV The total coupling strength, or 
lattice relaxation energy, is = 0.1 eV. Phonons above 
250 cm -1 give negligible contribution to As a conse- 
quence, they are not reported in the Table nor discussed 
here. As already mentioned in Section |IV[ the QHLD 
internuclear potential does not include the CT interac- 
tion, so that values in Table [T] are the zero-order refer- 
ence frequencies in the absence of Peierls coupling. In 
other words, they correspond to an hypothetical state 
with Ed = or, equivalently, to a state where the elec- 
tronic excitations are moved to infinite energy.^ As a 
consequence, the frequencies of Table [I] cannot be di- 
rectly compared to the experimental frequencies. In fact, 
phonons with high values of coupling constants are ex- 
pected to be the most intense in the spectra, and to occur 
at frequencies lower than the zero-order ones. Therefore 
we can for instance anticipate that the v§ and v% modes 
of Table [I] correspond to the group of intense bands oc- 



TABLE I: Calculated low-frequency B u phonons of TTF-CA 
in the N phase, cui, gi, and Si denote the zero-order frequency, 
Peierls coupling constant, and coupling strength, respectively. 



B u 



mode 


uo i (cm -1 ) 


gi (meV) 


Si (meV) 


v\ 


27.6 


-5.16 


7.8 


V2 


38.6 


8.30 


14.4 


^3 


55.8 


-8.68 


10.9 


Z/ 4 


80.5 


-2.16 


0.5 


V5 


90.5 


-9.66 


8.3 


V6 


99.4 


-16.99 


23.4 


V 7 


113.2 


-2.27 


0.4 


"8 


118.5 


-20.71 


29.2 


Vq 


134.0 


0.08 


0.0 


VlO 


194.3 


7.85 


2.6 


V\\ 


206.0 


-2.43 


0.2 


V\1 


211.3 


-2.93 


0.3 


V13 


250.7 


-13.75 


6.1 



curring around 70 cm -1 (Fig. |4|. Incidentally, the eigen- 
vectors of the vq phonon correspond to rigid molecular 
displacements along the stack, whereas the v% phonon 
has a more complex description, being a mixture of TTF 
torsional motion and molecular displacement along the 
stack. Modes above 200 cm -1 , on the other hand, are 
almost pure intramolecular modes, like for instance the 
^13 phonon, which essentially corresponds to a TTF out- 
of-plane motion. Their intensity may then have a signifi- 
cant "intrinsic" contribution (namely, not due to Peierls 
coupling), that is not accounted for in the present model 
calculation. 

Both experiment and calculations suggest that several 
modes are appreciably coupled to the CT integral, lead- 
ing to a multi-mode Peierls coupling. The resulting prob- 
lem becomes complex because the phonons are mixed 
through their common interaction with the CT electrons. 
Thus a meaningful comparison with experiment requires 
detailed modeling. The effect of Peierls coupling can be 
dealt with separately for the phonon frequencies and for 
their IR intensity.^ Accordingly, the analysis will be car- 
ried out in two steps. 

It has been shown that in the presence of Peierls 
electron-phonon interaction, the squared perturbed fre- 
quencies Q/j and normal modes coordinates Qj are ob- 
tained from the diagonalization of the following force 
constant matrix, written in the basis of the reference co- 
ordinates 



Fij — ujiUjjSij 



y/ViVj gigj Xb- 



(2) 



In this equation, uji and g^ are the reference frequencies 
and coupling constants of Table [l| and Xb is the electronic 
response to phonon perturbation. The Xb values needed 
to apply Eq. [2] are calculated as described in Ref. [13| by 
tuning the parameters of the modified Hubbard model as 
to mimic the behavior of TTF-CA. Since in going from 
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TABLE II: Calculated electronic response Xb and squared po- 
larizability derivative (dP/dS) for different degree of ionicity 
g. The correspondence between g and T is derived from ex- 
periment fRef. 120)) . 



Xt (eV- 1 ) 


(dP/35) 2 


Q 


T(K) 


5.7971 


0.0858 


0.197 


277 


6.1230 


0.1034 


0.209 


247 


6.4862 


0.1257 


0.222 


215 


6.8944 


0.1544 


0.236 


183 


7.3561 


0.1919 


0.250 


153 


7.8837 


0.2415 


0.266 


123 


8.4927 


0.3089 


0.282 


81 



200 to 82 K th e TTF-CA ionicity g changes from about 
0.2 to 0.3 j 23 ^ in the first column of Table || we list the 
Xb values corresponding to different ionicities in this in- 
terval. By matching the calculated ionicities in the third 
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we 



column of the Table with the experimental ones, 
can translate the ^-dependence of the calculated Xb m a 
T-dependence, as reported in the fourth column of Table 
\H\ In this way we can perform a more direct comparison 
with experimental data. 

Fig. [9] summarizes the temperature dependence of the 
perturbed frequencies Qj below 100 cm -1 calculated 
through Eq. with the parameters listed in Tables [i] 
and|ll{ The agreement between the calculated frequencies 
in Fig. [9] and the experimental ones presented in Fig. [6] 
is very good. First of all, the number of modes below 
100 cm -1 is the same, namely six. Second, an appre- 
ciable softening is detected only for the lowest frequency 
mode, and in the proximity of the phase transition. The 
weak softening shown by some higher frequency modes in 
Fig. [9] is likely compensated by the usual frequency hard- 
ening by lowering T, a factor not included in the model. 

The calculations quite naturally explain why - al- 
though several phonons are coupled to the CT electrons 
- the red shift of the associated bands on getting close 
to the transition is so small. As mentioned above, all 
the phonons are coupled together through their common 
interaction with the electronic system. Then, when two 
phonon frequencies get closer due to a softening of the 
higher frequency phonon, there is mixing in the phonon 
description, and the softening is "transferred" to the 
lower frequency phonon. In more precise terms, with low- 
ering temperature Xb increases with g (Table [il]), leading 
to an evolution of the normal modes Qj\ 



(3) 



where lij are the eigenvectors obtained, for each tem- 
perature, by diagonalizing the F matrix of Eq. [2] The 
Peierls coupling constants are accordingly modified: The 
coupling constants Gj in the basis of the perturbed nor- 
mal modes Qj are linear combination of the reference 
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FIG. 9: Temperature dependence of the Peierls-coupled 
modes, calculated by Eq. ([2]), on the basis of Table [I] and 
Table [IT] parameters. 



coupling constants: 




(4) 



and therefore evolve with Xb- As a result, on approach- 
ing the phase transition the coupling strength is progres- 
sively transferred to the lower frequency modes. Only in 
close proximity of the transition all the coupling strength 
collapses to the lowest frequency phonon, yielding sub- 
stantial softening of this specific mode. 

We now turn our attention to the intensity of the 
phonon bands. For a regular stack (equispaced D and 
A units along the stack, like N TTF-CA) the IR inten- 
sity of the Peierls modes is largely dominated by a term 
accounting for the charge fluctuations induced by oscilla- 
tion in the dimerization amplitude S. For a single Peierls 
mode with frequency cop the total oscillator strength isP^ 



/p 



m e d 2 LUped ( DP 
T 2 [~d5 



(5) 



where m e is the electronic mass, d the equilibrium dis- 
tance between D and A molecules, and P the electronic 
polarizability per site. The (dP/dS) 2 values calculated 
for parameters relevant to N TTF-CA are reported in 
the second column of Table [ll[ The temperature depen- 
dence of the (dP/d5) 2 terms compares well with the ex- 
perimentally determined growth of the spectral weight in 
the range of 0-100 cm -1 (Fig. 

In the multi-mode Peierls coupling case, the oscilla- 
tor strength is partitioned among the coupled modes. 
Specifically, each normal coordinate Qj modulates the 
CT integral as described by the coupling constants Gj in 
Eq. Q. Using the usual chain-rule, the (5-derivative of P 
m Eq. can be rewritten as a sum of Qj derivatives. 
After some algebra we derive the following expression for 



the IR oscillator strength of each mode coupled to the 
electronic degrees of freedom: 



fi 



m e d 2 2 ( dP 
t 2 j \ dS 



(6) 



We can now use the standard expression for the 
frequency-dependent dielectric constant 



6M 



3.96857-10 



10 N V 
V ^ 



(Si) 



- co 2 ) — iojTj 



(7) 



to calculate the contribution of the Peierls coupled modes 
to the low-frequency spectra of TTF-CA. The numerical 
factor in Eq. ^ refers to N/V (number of molecules per 
unit-cell volume) expressed in A -3 , and the frequencies 
oj, and damping Tj : in cm -1 . 

The calculated reflectivity is compared to the experi- 
mental one in Fig. [I] of Section III A The unperturbed 
frequencies and coupling constants needed for the calcu- 
lation are taken from Tab.[T| whereas Fj is set equal to 4.0 
cm -1 for all the modes. The values of Xb — 5.6667 eV _1 
and (dP/dS) 2 = 0.07922 are obtained from interpolation 
of the data in Tab. II In addition, we use d = 3.7 A, 
N/V = 4.92386 • lO^A" 3 (Ref. [26]), and t = 0.2 eV. 
The CT electronic transition is added in the calculation 
as an extra term in Eq. ([7]), where /ct, ^ct and Tct 
are derived directly from experiment.^ Finally, we put 
Coo = 2.5 to account for the high-energy contributions. 

Fig. [T] shows a very good agreement between model 
calculation and experiment, considering all the approx- 
imations involved in the model and in the estimate of 
the parameters, in particular of the gi. The frequen- 
cies are somewhat off, but well within the usual errors of 
QHLD calculations (± 10 cm -1 ). The calculated oscil- 
lators strengths for the lower frequency modes are just a 
little below the experiment. The agreement for the oscil- 
lator strength of the phonons at higher frequencies is not 
expected to be exact, because these phonons have intra- 
molecular character and, as discussed above, may have a 
non-negligible "intrinsic" intensity. The simulation indi- 
cates that at room temperature the most strongly cou- 
pled modes cluster around approximately 70 — 80 cm -1 . 



The discussion in Section [TTT] puts in evidence a quite 
complex evolution of TTF-CA spectra with temperature. 
Such behavior is qualitatively reproduced by the calcu- 
lation. An overall increase of reflectivity by lowering T 
is brought in by the increase in (dP/d5) 2 , and the shift 
of spectral weight towards lower frequency is accounted 
for by the change in the coupling constants Gj, as dis- 
cussed after Eq. Q. However, we can gain a better in- 
sight into the factors affecting the temperature evolution 
of the spectra if better agreement between experiment 
and simulation is achieved. To such aim, Tab. [I] frequen- 
cies and coupling constants are used as starting values 
for a nonlinear fitting procedure based on Eqs. Q and 
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FIG. 10: Experimental (thick line) and calculated (thin line) 
low- frequency reflectivity spectra of TTF-CA, parallel polar- 
ization, at three different temperatures. 
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FIG. 11: (color online) Calculated ratio of reflectance of the 
spectra below 110 K to the spectrum at 110 K. 



The variation of the frequencies and of the oscillator 
strengths with increasing g (lowering T) is induced in 
the model by the corresponding increase in both Xb an d 
(dP/dS) 2 , as detailed in Table [5] However, it is clear 
from experiment that also the bandwidths (damping) Tj 
change with T. In the model, the bandwidths are addi- 
tional parameters, but in the fit we choose to correlate 
them to the anharmonicity introduced by the electron- 
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phonon coupling. Following Fano, 2 ^ we impose the fol- 
lowing relationship to Tj (in cm -1 ): 



2.0 = k rjj + 2.0. 



(8) 



On top of the intrinsic linewidth of 2.0 cm -1 , the 
electron-phonon coupling gives an additional contribu- 
tion proportional to rjj = G^/flj. The proportionality 
factor k adds just one extra adjustable parameter to the 
fit, and Eq. ([8| represents a strong constrain. The fit 
shows a high sensitivity of spectra on the adjustable pa- 
rameters. Thus a variation of the unperturbed frequen- 
cies ou i or of the coupling constants gi by a few per cent, 
may profoundly alter the shape of the calculated spec- 
tra, especially in the case of the strongly coupled modes 
which cluster around 70 — 80 cm -1 . For this reason, it is a 
challenge to obtain a satisfactory fit at all temperatures. 



Fig. 10 illustrates the results of the fit of the reflectivity 
spectra at three temperatures. The fit is restricted to the 
spectral region below 150 cm -1 , so that the CT transition 
and its variation with temperature 2 ^ is simulated by . 
Furthermore, we allow the frequencies uoi to increase by 
a few wave number on lowering T, in order to simulate 
the usual hardening due to thermal contraction. The 
coupling constants and the proportionality factor for the 
damping, k, are instead kept temperature independent. 

Albeit not all details of the spectrum are precisely re- 
produced, the essential features are caught: By lowering 
the temperature the reflectivity increases, with a redistri- 
bution of the oscillator strengths and of the bandwidths 
as the spectral weight shifts towards lower frequencies. 
The downshift of spectral weight becomes more dramatic 
just before the phase transition. This is clearly demon- 
strated in Fig. [TT] which resembles the experimental be- 
havior reported in Fig. [3] The shift of the spectral weight 
is accompanied by an increase of the bandwidth of the 
^i, lowest frequency phonon, which becomes overdamped 
around 100 K and below (Fig. 10). The experimentally 
observed overdamping of the v\ phonon in Fig. [8] is then 
explained by the transfer of the entire electron-phonon 
coupling strength to this phonon just before the phase 
transition, with associated increase of anharmonicity. 



VI. DISCUSSION AND CONCLUSIONS 

In this paper we have presented temperature depen- 
dent polarized reflectivity spectra collected down to 
5 cm -1 for TTF-CA crystals above the NIT. A detailed 
theoretical analysis of the data has allowed us to offer the 
first direct identification of the Peierls softening mecha- 
nism. 

An indirect experimental evidence of the softening of 
lattice vibrations in TTF-CA was offered by the occur- 
rence of combination (two-phonon) bands in vibrational 
spectra collected in the region of the intramolecular vi- 
brations (mid-IR) ™ Several broad bands were found in 
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FIG. 12: (color online) Experimental (taken from Ref. I20[) 
and calculated frequency (squares) of the effective soft mode. 



IR spectra of TTF-CA that could not be assigned to fun- 
damental modes. These bands, whose intensity increases 
on approaching NIT, are symmetrically located at the 
low and high frequency side of a Raman band assigned to 
a totally-symmetric molecular vibration. The two side- 
bands are then assigned to the sum and difference combi- 
nations of the totally-symmetric molecular vibration with 
a lattice phonon. The frequency of this lattice phonon 
lowers from about 70 cm -1 to about 20 cm -1 when go- 
ing from room temperature to 81 K, as illustrated by 
In the light of the present work, in- 



the dots in fig. 12 



dicating a multi-mode Peierls coupling, the soft-phonon 
inferred from the analysis of combination bands must cor- 
respond to an "effective" Peierls phonon, resulting from 
the weighted contribution of the several Peierls-coupled 
modes. 

Specifically, the side-bands in the mid-IR region re- 
sult from the combination of a totally symmetric molec- 
ular vibration with a superposition of several B u modes. 
Thus the apparent peak frequency of the effective soft 
mode, fi e ff, is the weighted average of the frequencies of 
the Peierls phonons. In order to calculate fi e ff from the 
spectral simulation discussed in the previous Section, we 
adopt the following expression: 



^eff 



(9) 



where each frequency is weighted by rjj, the same factor 
entering the definition of vibrational linewidth in Eq. ([8| , 
and connected to the strength of the Peierls coupling. 
The computed value of ft e ^ is compared with experi- 
ment in Fig. [12] The very good agreement between the 
effective phonon frequency obtained from the analysis of 
combination bands in the mid-IR region and of lattice 
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phonons in the far-IR region strengthen our interpreta- 
tion of the spectroscopic data in terms of a soft mode. 
Quite interestingly, the same effective soft mode of TTF- 
CA quantitatively accounts for the peaks observed in the 
diffuse X-ray scattering dataP^ 

The softening of the lattice phonons, combined with 
the increase of their intensity on approaching NIT leads 
to very large vibrational contributions to the dielectric 
constant! 12 * 13 * By using the parameters found in the fit 
of reflectivity spectra, we can also estimate the tempera- 
ture evolution of the static dielectric constant, ei, which 
corresponds to the zero-frequency real part of the dielec- 
tric constant in Eq. [7| As shown in Fig. [l3j the present 
estimate of e\ agrees very well with available experimen- 
tal data, strongly supporting the vibrational origin of the 
dielectric anomaly at NIT, as due to the large charge os- 
cillations associated with the Peierls modes 

Data analogous to the present ones have been previ- 
ously reported Okimoto et a/.p^who discussed the reflec- 
tivity spectra along the stack of TTF-QBrCl 3 (QBrCl 3 : 
2-bromo-3,5,6-trichloro-j9-benzoquinone), a CT crystal 
similar to TTF-CA that undergoes the NIT with con- 
tinuous evolution of ionicity. However, in that work the 
reflectivity data have been collected only from 650 down 
to 25 cm -1 , so that the subsequent KKT transformation 
is delicate, and the discussion of conductivity spectra out- 
side the experimentally accessed region should be taken 
with caution. In any case a broad band was enucleated 
from conductivity spectra which softens, without major 
broadening, from about 60 cm -1 at 293 K to about 10 
cm -1 just above the critical temperature (71 K). This 
band was not ascribed to a soft mode, but rather to the 
pinned mode of the so-called neutral-ionic domain walls 
(NIDW)P Whereas more extensive measurements are 
required to fully address this issue, our work on TTF- 
CA sheds doubts on this interpretation. 

NIDWs were introduced theoretically by NagaossP as 
charged boundaries of / dimerized domains excited in the 
host N regular chain. Quite interestingly, they were orig- 
inally restricted to the close proximity of discontinuous 
NIT, while the energy of the domains is often too large 
to allow for thermal population. 29 However, NIDWs have 
been quite often invoked to explain the rich phenomenol- 
ogy of NIT In particular, the above mentioned dielec- 
tric anomaly and the combination bands in mid-IR spec- 
tra were initially attributed to NIDWs as well as the 
peak in the diffuse X-ray signal.^ As discussed above, 
this paper gives definitive support to a different picture 
that, without invoking exotic excitations, quantitatively 
explains the rich and variegated phenomenology of NIT 
to the increase of the effective Peierls coupling on ap- 
proaching NIT. In this way different and apparently un- 
related phenomena are all naturally and quantitatively 
explained in terms of the softening of lattice phonons 
coupled via a Peierls mechanism to delocalized electrons. 

The definitive and unambiguous support to the above 
picture is obtained by a careful analysis of the far-IR 
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FIG. 13: (color online) Experimental (taken from Ref. 1 10|) 
and calculated [Eq. |7|)] static dielectric constant of TTF-CA 
as a function of temperature. 



spectral region to directly identify the Peierls modes. The 
analysis is non-trivial: apart from the requirement for 
very refined experimental data extending down to very 
low- frequency, the spectral interpretation is based on a 
model for multi-mode Peierls coupling that combines lat- 
tice dynamics (QHLD) calculations with the modified 
Hubbard-model for the electronic structure. In fact, in 
the presence of several phonons the multi-mode Peierls 
softening is a complex phenomenon: the mode descrip- 
tion changes on approaching the transition, with associ- 
ated mode mixing, intensity redistribution and broaden- 
ing. In the course of this mixing, the softening "jumps" 
from one mode to the nearest one below, until in the prox- 
imity of the phase transition all the softening is trans- 
ferred to the lowest frequency mode. The mechanism, 
evident in the simulation presented in Fig. [9| explains 
the softening observed in Fig. [6] and the temperature 
behavior of mode intensities and linewidths in Fig. [7[ 
The increase of the electronic susceptibility as the sys- 
tem is driven towards the NIT implies an increased effec- 
tive coupling, and accordingly the spectral weight shifts 
towards zero frequency for T < 100 K. Indeed, in this 
temperature range the Peierls coupling is almost com- 
pletely transferred to the lowest-frequency mode. At the 
same time, the relevant bandwidth broadens, finally lead- 
ing to an overdamped behavior. For this reason, in the 
spectra we just see the reflectivity increase towards zero 
frequency (Figs. 2]and[3|, without being able to actually 
detect the full absorption. Overdamping has been often 
observed in ferroelectric phase transitions.^ In the case 
of TTF-CA, we are able to provide the microscopic origin 
of this effect. 
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